A GEOMETRIC MEASURE-TYPE REGULARITY CRITERION FOR 
SOLUTIONS TO THE 3D NAVIER-STOKES EQUATIONS 

Z. GRUJIC 

O 

Abstract. A local anisotropic geometric measure-type condition on the super- 
$-H ■ level sets of solutions to the 3D NSE preventing the formation of a finite-time sin- 

CL| gularity is presented; essentially, local one-dimensional sparseness of the regions of 

, intense fluid activity in a very weak sense. 
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1. Prologue 

Rigorous study of geometric depletion of the nonlinearity in the 3D Navier-Stokes 
equations (3D NSE) was initiated by Constantin in ||Co94| ; the approach was based 
on the singular integral representation formula for the stretching factor in the evo- 
lution of the vorticity magnitude featuring a geometric kernel depleted by coher- 
ence of the vorticity direction. This representation was subsequently utilized by 
Constantin and Fefferman in IICoFe93 1 to show that as long as the vorticity direc- 
tion is Lipschitz-coherent, no finite-time blow up can occur and later by Beirao 



da Veiga and Berselli in | daVeigaBe02 1, where the Lipschitz-coherence regularity 
condition was replaced with ^ -Holder. 

Spatiotemporal localization of the i -Holder-coherence regularity criterion was 
performed in |GrZh06, Gr09|, and also - utilizing a different localization proce- 
dure - by Chae, Kang and Lee in IIChKaLe07l . 

The i -Holder-coherence condition is super-critical with respect to the natural 
scaling of the 3D NSE; a family of scaling-invariant, critical, hybrid geometric- 
analytic local regularity criteria - including a scaling-invariant improvement of 
the i -Holder-coherence condition - was presented in IGrGulO-lll . 

In the realm of the mathematical theory of turbulence, the i -Holder-coherence 
condition was recently (IDa Grll-lJ ) paired with the condition on a modified Kraich- 
nan scale to obtain a first rigorous evidence of existence of (anisotropic) enstrophy 
cascade in 3D viscous incompressible flows. 

A different approach to discovering geometric scenarios ruling out formation 
of singularities in the 3D NSE was introduced in [GrOlJ. The main idea was to 
utilize the local-in-time spatial analyticity properties of solutions in U' ( ||GrKu98l ) 
via the plurisubharmonic measure maximum principle - a generalization to C" (cf . 
llSadSll ) of the classical harmonic measure majorization principle in the complex 
plane (the log-convexity of the modulus of an analytic function; see, e.g., IINev70| ). 

The regularity criterion derived in BGrOll is a condition on the regions of in- 
tense fluid activity near a possible blow up time requiring local existence of a 
sparse coordinate projection on the scale comparable to the uniform radius of spa- 
tial analyticity. The estimate on the plurisubharmonic measure was performed 
within the framework of product-type domains - hence the requirement on a co- 
ordinate projection. This could be somewhat relaxed, but not substantially due 
to the rigidity of the C" structure. Also, once the computation of the plurisub- 
harmonic measure was reduced to the computation of the individual (coordinate) 
harmonic measures, the estimate on the harmonic measure was carried out with 
respect to an infinite strip, giving the argtmient a nonlocal character. 

In the present work, we completely bypass the rigidity of the C" structure, re- 
sulting in a much weaker local geometric measure-type condition. Utilizing trans- 
lational and rotational invariance of the 3D NSE, as well as some basic geomet- 
ric properties of the harmonic measure, the argument is ultimately reduced to 
the problem of estimating the harmonic measure of an arbitrary closed subset of 
[— 1, 1] computed at with respect to the unit disk. This is a generalization of 
the classical Beurling's problem IIBeu33[ [NevZOI proposed by Segawa in |Seg88|; 
a symmetric version was solved by Essen and Haliste in [EssHa89 1 , and the gen- 
eral case relatively recently by Solynin in [Sol99J via a general symmetrization 
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argument. The utility of Solynin's result here is that it allows us to formulate the 
condition in view in terms of the ratio of one-dimensional Lebesgue measures, 
rather than in terms of a specific 'rearrangement' of the one-dimensional trace of 
the region of intense fluid activity. The engine behind the proof is the interplay be- 
tween the diffusion in the model - quantified by the local-in-time (sharp) analytic 
smoothing in L°° - and the geometric properties of the harmonic measure. 

A precise statement of our regularity criterion will be given in the main text; at 
this point, we convey the essence of the result. Denoting the region of intense fluid 
activity at time s by Q,s{M) and the lower bound on the uniform radius of spatial 
analyticity by p{s) (this take places near a potential singular time), the condition 
in view is simply a stipulation that for a given point a;o/ there exists a radius r — 
r{xo), < r < p{s) and a unit vector d = d{xo), such that 

(1-1) < ^ 

for some 5 in (0, 1). (It will transpire that it is enough to require dl.ll l on a suitably 
chosen finite sequence of times .) There are two versions of the result, one for 
the velocity and one for the vorticity, based on the spatial analyticity estimates on 
solutions to the velocity and vorticity formulations of the 3D NSE, respectively. 

It is plain that 111. lb is a much weaker condition than the one in [GrOl]; all that is 
needed here is local sparseness of a one-dimensional trace of the region of intense 
fluid activity in a very weak sense. On the other hand, it is of a different nature 
and hence not directly comparable to the coherence of the vorticity direction-type 
regularity criteria. 

The main efficacy of the aforementioned regularity criterion is in ruling out 
various sparse geometric scenarios for a finite-time blow up, both in the velocity 
and the vorticity formulations. As succinctly put by P. Constantin, "intermittency 
implies regularity" |PC11|. An example of interest that can be ruled out is a blow 
up scenario in which the region of intense vorticity (defined as a region in which 
the vorticity magnitude exceeds a suitable fraction of the i°°-norm) is - at suitable 
near-blow up times - comprised of vortex filaments with diameters of the cross- 
sections boimded above by the uniform radius of spatial analyticity. 

The paper is organized as follows. In Section 2, we collect relevant properties 
of the harmonic measure in the plane, and in Section 3 we recall the local-in-time 
spatial analyticity of solutions in L°°. Section 4 contains the main result. The last 
section indicates a scenario which - in a statistically significant sense - leads to 
closing the scaling gap in the regularity problem. 

2. Harmonic measure 

Basic properties of the harmonic measure in the complex plane can be found, 
e.g., in liNevZOl lAhllOI . First, we briefly recall a few relevant facts following 
IIAhllOI . 

Let 57 and K be an open and a closed set in the complex plane, respectively. 
When the geometry of $7 \ if is not too convoluted, there exists a unique bounded 
harmonic function on H \ K, denoted by w = fl, K), such that - in the sense 
of a well-defined limit as a point approaches the boundary - w is equal to 1 on if 
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and on the rest of the boundary; ijj{z, il, K) is the harmonic measure of K with 
respect to VL computed at z. 

Two straightforward consequences of the general harmonic measure majoriza- 
tion principle, c.f. Theorem 3.1 MAhllOl , are the following ( MAhllOl , p. 39). 

Proposition 2.1. The harmonic measure is increasing (as a measure) with respect to both 
K and Q. 

Proposition 2.2. Let f be analytic in fl \ K, \ f\ < M, and \ f\ < m on K (in the sense 
of lim sup as a point approaches the boundary). Then 

\f{z)\ < m'^M^-'^ 

for any zinVL \ K, where 6 = uj{z, fl, K). 

This a refined form of the maximum modulus principle for analytic fimctions 
in \ A' (the log-convexity of the modulus of / - sometimes referred to as "two- 
constants theorem"). 

Another useful property of the harmonic measure is the following (see., eg., 
LNev70.n . 

Proposition 2.3. The harmonic measure is invariant with respect to conformal mappings. 

Finally, we recall a result on extremal properties of the harmonic measure in the 
unit disk B obtained by Solynin in | Sol99 1 . 

Theorem 2.1. Let K be a closed subset of [—1, 1] such that \K\ — 2 A /or some A, < 
A < 1, and suppose that € D \ AT. Then 

a;(0, D, K) > w(0, D, Kx) = - arcsin | ' j| ~ ^j' 

TT 1 + (1 — A)"^ 

where Kx = [-1, -1 + A] U [1 - A, 1]. 

The above theorem provides a generalization of the classical Beurling's result 
l|Beu33 l in which K is a finite union of intervals lying on one side of the origin. This 
was conjectured by Segawa in | Seg88J , and the symmetric version was previously 
resolved in |EssHa89| . 



3. Spatial analyticity in L°° 
The 3D NSE equations read 
(3.1) Mf + (it • V)u = -Vp + Am 

supplemented with the incompressibility condition Wu = 0, where u is the velocity 
of the fluid and p the pressure (the viscosity is set to 1). 

A method for deriving explicit local-in-time lower bounds on the uniform ra- 
dius of spatial analyticity of solutions to the NSE in Lp was introduced in |GrKu98l ; 
see also l|Ku991 for analogous results in the vorticity formulation. We will make 
use of the following sharp analyticity estimate in L°° (cf. BGulOl ; l|Ku03] for the 
corresponding real result). 
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Theorem 3.1. Let uq be in L°°{M?). Then, there exists an absolute constant co > 1 such 
that setting T = — r-^r-, a unique mild solution u = u(t) on [0, T] has the analytic 

4\K\\lo 

extension U = U{t) to the region 

TZt = {x + iy eC^ : \y\ < —Vt} 

Co 

for any t in (0, T]. In addition, 

\\U{t)\\L'^(n,) < collwolloo 

for all tin [0,T]. 

Recall that the vorticity formulation of the 3D NSE reads 

(3.2) + (u • V)w = (w • V)u + Auj 

where co ~ curl u is the vorticity. 

The vorticity version of the above theorem is as follows (the proof is analogous; 
utilizing the Biot-Savart law to close each iteration). 

Theorem 3.2. Let ojq be in (R^). Then, there exists an absolute constant do > I such 
that setting T = — - — , a unique mild solution u = uj{t) on [0, T] has the analytic 

"oll'^olloo 

extension Q, = Ct{t) to the region 

TZt = {x + iyeC^ : \y\ < ^Vi} 

do 

for any t in (0, T]. In addition, 

\\^{t)\\L-^{nt) < c!o||'^o||oo 

for all tin [0,r]. 



4. The main result 

We start with introducing a geometric measure-theoretic concept of weak local 
linear sparseness of a set around a point, at a given scale, suitable for our purposes. 



Definition 4.1. Let xq be a point in R^, r > 0, S an open subset ofR^ and 6 in (0, 1). 

arse around xq at scale 

15 n (a;o — rd, xo + rd)\ 



The set S is linearly J-sparse around xq at scale r in weak sense if there exists a 
unit vector d in S'^ such that 



2r 



< 5. 



In what follows, we derive the main result for the velocity formulation and sim- 
ply state the analogous result for the vorticity formulation; modifying the proof in 
the second case is essentially relabeling. 

For M > 0, denote by Ot(M) the super-level set at time t; more precisely, 

Ot(M) = {a; e : \u{x,t)\ > M}. 

Then, our main result reads as follows. 
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Theorem 4.1. Suppose that a solution u is regular on an interval (0, T*). (Recall that u 
is then necessarily in C^(0, T*); L°°^ .) 

Fix 5 in (0, 1), and let h — h{5) = ^ arcsin and a — a{5) > ^-j^. Assume that 
there exists e > such that for any t in (T* — e, T), either 

(i) t + on ! Mio > T* (cq is the constant featured in Theorem \5l\ , or 

(ii) there exists s — s{t) in t + 4c^||J(t)||2 + ^-Jn^^^^p such that for any spatial 
point xq, there exists a scale r,0 < r < 2c'^\]u{t)\\ — ' ^'^^ property that the super-level 
set ris(M) is linearly S-sparse around xo at scale r in weak sense; here, M = M{5) = 
^l|w(i)l!oo. 

Then, there exists 7 > such that u is in L°° {{T* - e, T* + 7) ; L°°) , i.e., T* is not a 
singular time. 

Proof. There are two cases to consider. 

Case I (i) holds for some t in (T* — e, T*). 

\n this case, the statement of the theorem follows from Theorem |3.1[ setting the 
initial time to t. 

Case II (ii) holds for all t in (T* - e, T*). 

Pick a time io in {T* — e,T*), and let S,h,a,M and sq = s(io) be as in the 
statement of the theorem. Then, for any xq in M^, there exists an r = r{xo) , < 
T < , 2|| }f ^„ — , and a direction vector d ~ d(xo), such that 

.sCq II ti^^to j II oc 

\n,„{M)nixo-rd,Xo+rd)\ ^ ^ 
2r 

The intent is to show 

(4.1) ll"(so)||oo < \Hto)\\oc. 

Fix xq. Recall that the NSE exhibit translational and rotational invariance in 
the spatial variable. Translate for —xq, rotate by the matrix Q transforming the 
unit direction d to the coordinate vector ei and denote the transformed solution 
by Ux„^Q. Then, Uxo,Qix,t) = Qu[Q^^{x + XQ),t). 

Solve the NSE locally-in-time starting at to ; the spatial analyticity properties of 
the solution at time Sq are given by simply translating in time the statement of 
Theorem |3.1l 

Moreover, since the rotation Q has no effect on computing the norms, the trans- 
formed solution Uj;g_Q at time sq enjoys exactly the same analyticity features. 

In particular - focusing on the first coordinate - Uxg^qiso) is spatially analytic 
on a strip, symmetric around the real axis, with the width equal to (at least) 

^(■^o) " 211 u Ml ■ 

The region of interest is the disk around the origin with the radius r, D^. Note that 
Dr is contained in the domain of analyticity of Uxo,q{so). 
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Our goal is to obtain an improved estimate on Uxo,q{0, sq). Denote by K the 
complement of the image of the set flsf, (M) n (xo — rd, xq + rd), under the change 
of coordinates, in [— r, r] . Then, K is closed, and the sparseness assumption implies 
\K\ > 2r(l - 6). li e K, \uxg^Q{0, sq)\ < \\u{to)\\oo, and we are done (with this 
a;o). If not, the harmonic measure maximum principle - Proposition l2.2l - together 
with the _L°°-bound on the complexified solution stated in Theorem 13. 1[ implies 



(4.2) |w,„,q(0,So)| < (— ||u(to)||ooj [co\Hto)\\oo) 

Recall that the harmonic measure is invariant under conf ormal mappings (Propo- 
sition |Z3] l. In particular, it is invariant under the scaling map z i-> iz. This paired 
with the monotonicity of the harmonic measure with respect to K (Proposition 
12.1b and Theorem 12 . II yields 



2 1 - (5^ 

(4.3) a;(0, Dr, K)>- arcsin = h. 

TT 1 + 0^ 

Combining the estimates (|4.2| | and | |4.3b leads to 



(4.4) |u,„,q(0,so)| < (^ll"(io)||oo)''(co||ii(io)||oo)' < \Hto)\\^. 

This, in turn, implies |u(a;o, so)| < |lw(<o)|loo/ and since xa was an arbitrary spa- 
tial point in R^, ||u(so)||oo < ||u(io)||oo- 

Let 

Mo = \Hto)\\oo- 

Setting ii = sq and repeating the argument yields |iu(si)||oo < ||u(ii)||oo < Mq, 
where si = s{ti). After finitely many steps, we reach the time s„, s„ < T* such 

that ||u(s„)||oo < Mq and s„ H — ^ > T*. The statement of the theorem now fol- 

cgA/o 

lows from Theorem |3.11 setting the initial time to s„. □ 



Remark 4.1. It is plain from the proof that it is enough to assume the condition on 
a finitely many suitably chosen times. 

The vorticity version and the proof are completely analogous - utilizing Theo- 
rem |3.2| in place of Theorem 13. II 

For M > 0, denote by Qf{M) the vorticity super-level set at time t; more pre- 
cisely, 

n'^{M) = {a; e R3 : \u;{x,t)\ > M}. 
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Theorem 4.2. Suppose that a solution u is regular on an interval (0, T*). 

Fix 5 in (0, 1), and let h ~ h{S) — |- arcsin and a = a{5) > Assume that 
there exists e > such that for any t in (T* — e, T), either 

(i) t + ,„,, ^, ... — > T* (dn is the constant featured in Theorem \3^ , or 

d^\\uj{t)\\r^ 

(ii) there exists s — s[t) in t + 4^ij|| J^^^n — , t + — such that for any spatial 
point xq, there exists a scale r, < r < — ^, with the property that the super-level 

2rfgll<^(t)l|l= 

set VL'^{M) is linearly S-sparse around xq at scale r in weak sense; here, M = M{5) = 



Then, there exists 7 > such that uj is in L°°{{T* - e,T* + •y); 1°°^ i.e., T* is 
a singular time. 



not 



5. Epilogue 

Direct numerical simulations of turbulent flows reveal (see, e.g., [SJ091| ) that 
the preferred geometric signature of the regions of intense vorticity is the one of 
vortex filaments. The general agreement seems to be that the length of a filament 
is - in a statistically significant sense - comparable with the macro scale, while 
the scaling of the diameter of the filament's cross-section seems to be harder to 
pin down (although, it is mostly found to be comparable to some version of Kol- 
mogorov dissipation scale). For rigorous mathematical results concerning creation 
and dynamics of vortex tubes in turbulent flows, the reader is referred to ICPS95I . 

Let us for a moment adopt the aforementioned geometry as a blow up scenario, 
and define the region of intense vorticity at a near-blow up time t to be the region 
in which the vorticity magnitude exceeds a fraction of 1 1 w ) 1 1 00 • Then, Theorem l4.2l 
implies that as long as the diameters of the filaments' cross-sections are dominated 

by —, for a suitable constant Ci > 1, no blow up can occur. At this point, 

^1 Mt)\% 

recall that starting with the initial vorticity a finite Radon measure, the L^-norm 
of the vorticity is bounded - uniformly in time - over any interval (0, T)) |,Co90| . 
Tchebyshev then implies the decrease of the distribution fimction of the vorticity 
of at least i; consequently, the volume of the region of intense vorticity decreases 

at least as C2-r. — -p-r. — . Assuming that the length of a filament is comparable with 
the macro scale, this implies the decrease of the diameters of the filaments' cross- 
sections of at least C3 —, which is exactly the scale needed for the application 

^ Mt)\\lo 

of Theorem l4.2l 

The above ruminations offer a geometric scenario leading to closing the scaling 
gap in the regularity problem. Assuming that the 'shape', i.e., the general geom- 
etry is correct, the weakest link is the assumption that the length of a filament be 
comparable to the macro scale; this was simply borrowed from the picture painted 
by the numerical simulations. However, in a very recent work | DaGrll -2|, the au- 
thors utilized the ensemble averaging process introduced in their study of turbulent 
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cascades in physical scales of 3D incompressible flows ( IIDaGrlOII ) to arrive at pre- 
liminary results indicating that the averaged vortex-stretchrng term is - near a pos- 
sible blow up time T* - positive across a range of scales extending from a power 
of a modified Kraichnan scale to the macro scale. This provides a mathematical 
evidence of persistence of the macro scale-long vortex filaments (in a statistically 
significant sense), and the pertaining research will be pursued in the future. 
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